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We study a quantum point contact in a quantum spin Hall insulator. It has recently been shown 
that the Luttinger liquid theory of such a structure maps to the theory of a weak link in a Luttinger 
liquid with spin with Luttinger liquid parameters g p = 1/ g a = g < 1. We show that for weak 
interactions, 1/2 < g < 1, the pinch-off of the point contact as a function of gate voltage is 
controlled by a novel quantum critical point, which is a realization of a nontrivial intermediate fixed 
point found previously in the Luttinger liquid model with spin. We predict that the dependence of 
the conductance on gate voltage near the pinch-off transition for different temperatures collapses 
onto a universal curve described by a crossover scaling function associated with that fixed point. 
We compute the conductance and critical exponents of the critical point as well as the universal 
scaling function in solvable limits, which include g = 1 — e, g = 1/2 + e and g = l/\/3. These 
results, along with a general scaling analysis provide an overall picture of the critical behavior as a 
function of g. In addition, we analyze the structure of the four terminal conductance of the point 
contact in the weak tunneling and weak backscattering limits. We find that different components 
of the conductance can have different temperature dependence. In particular, we identify a skew 
conductance Gxy , which we predict vanishes as T 1 with 7 > 2. This behavior is a direct consequence 
of the unique edge state structure of the quantum spin Hall insulator. Finally, we show that for 
strong interactions g < 1/2 the presence of spin non conserving spin orbit interactions leads to a 
novel time reversal symmetry breaking insulating phase. In this phase, the transport is carried by 
spinless chargons and chargeless spinons. These lead to nontrivial correlations in the low frequency 
shot noise. Implications for experiments on HgCdTe quantum well structures will be discussed. 

PACS numbers: 71.10.Pm, 72.15.Nj, 85.75.-d 



I. INTRODUCTION 

A quantum spin Hall insulator (QSHI) is a time re- 
versal invariant two dimensional electronic phase which 
has a bulk energy gap generated by the spin orbit 
interactioni 2 -. It has a topological order— which requires 
the presence of gapless edge states similar to those that 
occur in the integer quantum Hall effect. In the sim- 
plest version, the QSHI can be understood as two time 
reversed copies of the integer quantum Hall stated for up 
and down spins. The edge states, which propagate in 
opposite directions for the two spins, form a unique one 
dimensional system in which elastic backscattering is for- 
bidden by time reversal symmetry 1 . This state occurs in 
HgCdTe quantum well structures 5 , and experiments have 
verified the basic features of the edge states, including the 
Landauer conductance 6 - 2e 2 /h, as well as the non locality 
of the edge state transport 7 -. 

In the presence of electron interactions, the edge states 
form a Luttinger liquid 8 - ' 9 ' 10 ' 11 ' 12 ' 13 ' 14 . For strong inter- 
actions (when the Luttinger liquid parameter g < 3/8) 
random two particle backscattering processes destabilize 
the edge states, leading to an Anderson localized phase. 
For g > 3/8 (or a sufficiently clean system), however, one 
expects the characteristic power law behavior for tunnel- 
ing of a Luttinger liquid. 

A powerful tool for probing edge state transport ex- 
perimentally is to make a quantum point contact. As 
depicted in Fig. QJa,b), a gate voltage controls the cou- 
pling between edge states on either side of a Hall bar as 
the point contact is pinched off. Recently, the point con- 



tact problem for a QSHI has been studie d 10 ' 11 . Hou, Kim 
and Chamoni - made the interesting observation that the 
QSHI problem maps to an earlier studied mode l 15 ' 16 of 
a weak link in a spinful Luttinger liquid (SLL), in which 
the charge and spin Luttinger parameters are given by 
g p = g and g a = l/g~. For sufficiently strong inter- 
actions (g < 1/2) they found that the simple perfectly 
transmitting and perfectly reflecting phases are both un- 
stable. They showed that as long as spin is conserved 
at the junction the low energy behavior is dominated by 
a non trivial "mixed" fixed point of the SLL, in which 
charge is reflected but spin is perfectly transmitted. This 
charge insulator/spin conductor (IC) phase leads to a 
novel structure in the four terminal conductance of the 
point contact. 

In this paper, we will focus on the QSHI point con- 
tact for weaker interactions, when 1/2 < g < 1. In this 
regime the open limit (or weak backscattering, "small v" ) 
and the pinched off limit (or weak tunneling, "small t v ) 
are both stable perturbatively. This is different from the 
behavior in an ordinary Luttinger liqui d 15 ' 16 ' 18 or a frac- 
tional quantum Hall point contac t 19 ' 20 . In those cases 
the perfectly transmitting limit is unstable for g < 1. 
Weak backscattering is relevant and grows at low energy, 
leading to a crossover to the stable perfectly reflecting 
fixed point. The fact that both the small v and the small 
t limits are stable for the QSHI point contact means that 
there must be an intermediate unstable fixed point which 
separates the flows to the two limits. This unstable fixed 
point describes a quantum critical point where the point 
contact switches on as a function of the pinch-off gate 
voltage. We will argue that in the limit of zero tempera- 
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FIG. 1: A quantum point contact in a QSHI, controlled by 
a gate voltage Vg- In (a) Vg < Vq, and the point contact 
is pinched off. The spin filtered edge states are perfectly 
reflected. In (b) Vg > Vq, and the point contact is open. 
The edge states are perfectly transmitted. In (c) we plot the 
conductance (later defined as Gxx) as a function of Vg for 
different temperatures. As the temperature is lowered, the 
pinch-off curve sharpens up with a width T a . The curves 
cross at a critical conductance G* , and the shape of the curve 
has the universal scaling form Ql.ip . The plotted curves are 
based on Eq. 13.91 valid for g = 1 — e, which is computed in 
Section III.C. 

ture the point contact switches on abruptly as a function 
of gate voltage Vg, with conductance G = for Vg < Vq 
and G = 2e 2 /h for Vg > Vq. At finite but low tem- 
perature T, the shape of the pinch-off curve G(Vg,T) is 
controlled by the crossover between the unstable and sta- 
ble fixed points, and is described by a universal crossover 
scaling function, 

lim G(V G ,T) = ^-g g (c^-). (1.1) 

Here AVg = Vg — Vq and c is a non universal constant. 
a g is a critical exponent describing the unstable inter- 
mediate fixed point. Q g (X) is a universal function which 
crosses over between and 1 as a function of X. a g and 
Q g {X) are completely determined by the Luttinger liquid 
parameter g. This behavior means that as temperature is 
lowered, the pinch-off curve as a function of Vg sharpens 
up with a characteristic width which vanishes as T Qs , as 
shown schematically in Fig. QJc). The curves at different 
low temperatures cross at G* = g {0), the conductance 
of the critical point. Eq. 11.11 predicts that data from 
different temperatures can be rescaled to lie on the same 
universal curve. 

The crossover scaling function Q g (X) is similar to the 
scaling function that controls the lineshape of resonances 



in a Luttinger liquid^ 6 ' 21 and in a fractional quantum Hall 
point contact^. That scaling function was computed ex- 
actly for all g by Fendley, Ludwig and Saleur— using the 
thermodynamic Bethe ansatz. That problem, however, 
was simpler than ours because the critical point occurs 
at the weak backscattering limit, which is described by 
a boundary conformal field theory with a trivial bound- 
ary condition 23 . The intermediate fixed point relevant to 
our problem has no such simple description. Thus, even 
the critical point properties a g and G* (which were sim- 
ple for the resonance problem) are highly nontrivial to 
determine. 

Intermediate fixed points in Luttinger liquid problems 
were first discussed in Refs. UH and [la in the context of 
SLLs. However, for that problem they occur in a rather 
unphysical region of parameter space g a > 2, because 
spin rotational invariance requires g a = 1. To our knowl- 
edge, the QSHI point contact provides the first physically 
viable system to directly probe these non trivial fixed 
points. 

The existence of the intermediate fixed points can be 
inferred from the stability of the simple perfectly trans- 
mitting or reflecting fixed point s 15 ' 16 . However their 
properties are difficult to compute, and a general char- 
acterization of these critical points remains an unsolved 
problem in conformal field theory 2 ^. Two approaches 
have been used to study their properties. In Ref. [H, 
a perturbative approach was introduced which applies 
when the Luttinger parameters are close to their criti- 
cal values g* and g*, where the simple fixed points be- 
come unstable. (For instance, for the weak backscatter- 
ing limit, g* = 1/2, g* = 3/2). For 5pj(T = g* pa ~ e, the 
fixed point is accessible in perturbation theory about the 
simple fixed point, and it's properties can be computed 
in a manner analogous to the e expansion in statistical 
mechanics. 

An alternative approach is to map the theory for spe- 
cific values of g p and g a onto solvable models. In Ref. 
l25l Yi and Kane recast the Luttinger liquid barrier prob- 
lem as a problem of quantum Brownian motion (QBM) 
in a two dimensional periodic potential. When g p — 1/3, 
g„ = 1 and the potential has minima with a honeycomb 
lattice symmetry, a stable intermediate fixed point which 
occurs in that problem was identified with that of the 
3 channel Kondo problem. This, in turn is related to 
the solvable SU(2) 3 Wess Zumino Witten model^, al- 
lowing for a complete characterization of the fixed point. 
This idea was further developed by Affleck, Oshikawa and 
Saleur—, who provided a more general characterization of 
the fixed point in terms of the boundary conformal field 
theory of the three state Pott's model. For g p = l/v3 
and g a — v3 the QBM model with triangular lattice sym- 
metry has an unstable intermediate fixed point, which we 
will see is related to the fixed point of the QSHI prob- 
lem. In Ref. [25] symmetry arguments were exploited to 
determine the critical conductance G* in that case. 

In this paper we will compute a g and Q g (X) (along 
with a multiterminal generalization of the conductance) 
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in three solvable limits: 

(i) For g = 1 — e, we will perform an expansion for 
weak electron interactions. For non interacting electrons 
the point contact can be characterized in terms of a scat- 
tering matrix S^j, for arbitrary transmission. Weak in- 
teractions lead to a logarithmic renormalization of Sij. 
Following the method developed by Matveev,Yue and 
Glazman 27 , this allows G g (X) and a g to be calculated 
exactly in the limit g — > 1. 

(ii) For g — 1/2 + e we find that the intermediate 
fixed point approaches the charge insulator/spin conduc- 
tor fixed point, allowing for a perturbative calculation of 
the fixed point properties G* and a g to leading order in 
e. Moreover, for g = 1/2 the Luttingcr liquid theory can 
be fermionized, which allows the full crossover function 
Qg(X) to be determined in that limit. 

(Hi) For g = l/VS the self duality argument developed 
in Ref. [25| allows us to compute the fixed point conduc- 
tance G* exactly. 

These three results, along with the general scaling 
analysis provide an overall picture of the critical behavior 
of the QSHI point contact as a function of g. 

In addition to the analysis of the pinch-off transition 
discussed above, we will touch on two other issues in this 
paper. First, we will introduce a convenient parameteri- 
zation of the four terminal conductance as a 3 x 3 conduc- 
tance matrix. In this form symmetry constraints on the 
conductance are reflected in a natural way. Moreover, we 
will predict that different components of the conductance 
matrix have different temperature dependence at the low 
temperature fixed points. In particular, we will introduce 
a "skew" conductance Gxy, which is predicted to vanish 
as T 7 with 7 > 2. For non interacting electrons we will 
show that Gxy = 0, and for weak interactions 7 = 2. 
This behavior is a direct consequence of the spin filtered 
nature of the edge states, and does not occur in a generic 
four terminal conductance device. It is thus a powerful 
diagnostic for the edge states. 

Secondly, we will examine the role of spin orbit terms 
at the point contact which respect time reversal symme- 
try but violate spin conservation. For g > 1/2 we will 
provide evidence that such terms are irrelevant at the in- 
termediate critical fixed point, so that they are unimpor- 
tant for the critical behavior of the point contact. How- 
ever, for g < 1/2, such terms are relevant. Hou, Kim and 
Chamon 10 pointed out that these terms are relevant per- 
turbations at the charge insulator/spin conductor fixed 
point for g < 1/2, but they did not identify the stable 
phase to which the system flows at low energy. We will 
argue that the system flows to a time reversal symme- 
try breaking insulating state in which the four terminal 
conductance = 0. Since spin orbit interaction terms 
will generically be present in a point contact, the true 
low energy behavior of a point contact will be described 
by this phase. An interesting consequence of the broken 
time reversal symmetry of this phase is that the weak 
tunneling processes which dominate the conductance at 
low, but finite temperature are not electron tunneling 



processes. Rather, they involve the tunneling of neutral 
spinons and spinless chargons. This has nontrivial impli- 
cations for four terminal noise correlation measurements. 
A related effect has been predicted by Maciejko et alJ^ 
for the insulating state of a single impurity on a single 
edge of a QSHI. This insulating state, however, requires 
stronger electron electron interactions. It occurs in the 
regime g < 1/4, where weak disorder already leads to 
Anderson localization. 

This paper is organized as follows. In section II we 
discuss our model and analyze five stable phases. In ad- 
dition to the simple fixed points, where charge and spin 
are either perfectly reflected or perfectly transmitted, we 
discuss the time reversal symmetry breaking insulating 
phase which occurs for strong interactions with spin or- 
bit. In section III we discuss the critical behavior of the 
conductance at the pinch-off transition. We will begin 
in section III. A with a general discussion of the scal- 
ing theory and phase diagram, along with a summary 
of our results. Readers who are not interested in the 
detailed calculations can go directly to this subsection. 
In the following subsections we describe the calculations 
for g — 1/V3, 9 = 1 — e and g = l/2 + ein detail. In 
section IV we conclude with a discussion of experimental 
and theoretical issues raised by this work. In appendix 
A we show describe our parameterization of the four ter- 
minal conductance and show that in this representation 
symmetry constraints have a simple form. 



II. MODEL AND STABLE PHASES 

In this section we will describe the Luttinger liquid the- 
ory of the QSHI point contact. We will begin in section 
II. A by describing the Luttinger liquid model first for a 
single edge and then relating the four edges to the the- 
ory of the SLL. We then discuss the four terminal conduc- 
tance. In section II. B we describe the simple limits of our 
model which correspond to stable phases. The simplest 
limits are the perfect transmission limit, or charge con- 
ductor/spin conductor (CC), the perfect reflection limit, 
or charge insulator/spin insulator (II). In addition we will 
discuss the "mixed" phases, including the charge insula- 
tor/spin conductor (IC) and the charge conductor/spin 
insulator (CI). 

For most of this section we will assume that spin is con- 
served. While spin nonconserving spin orbit interactions 
are allowed and will generically be present we will argue 
that they are irrelevant for the fixed points and crossovers 
of physical interest. An exception to this, however, oc- 
curs for strong interactions when g < 1/2. This will be 
discussed in section II. B. 5, where we will show that there 
are relevant spin orbit terms which destabilize the CC, 
II and IC phases. We will argue that these perturba- 
tions flow to a different low temperature phase, which 
we identify as a time reversal symmetry breaking insula- 
tor (TBI). In that section we will explore the transport 
properties of that state. 
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Much of the theory presented in this section is con- 
tained either explicitly or implicitly in the work Hou, Kim 
and Chamon^, as well as in Refs. Illlllalla . We include it 
here to establish our notation and to make our discussion 
self contained. We will highlight, however, three results 
of this section which are original to this work. They in- 
clude (1) our analysis of the four terminal conductance, 
which predicts that different components of the conduc- 
tance matrix have different temperature dependence. In 
particular, we find that the skew conductance Gxy van- 
ishes at low temperature as T 7 with 7 > 2. (2) In sec- 
tion II. B. 5 we introduce the TBI phase discussed above. 
(3) We introduce a perturbative analysis of the IC and 
CI phases in section II. B. 3 and II. B. 4. While this was 
partially discussed in Ref. HH, we will show that a full 
analysis requires the introduction of a pseudo-spin degree 
of freedom in the perturbation theory. This new pseudo- 
spin does not affect the lowest order stability analysis of 
the IC phase, but it will prove crucial for the second or- 
der renormalization group flows, which will be used in 
the e expansion in Section III.D. 



A. Model 

The edge states on the four edges in Fig. [lja,b) ema- 
nating from the point contact may be described by the 
Hamiltonian 
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Here V'i.in and ipi t out are a time reversed pair of fermion 
operators with opposite spin which propagate toward and 
away from the junction, vq is the bare Fermi velocity, 
and u is electron interaction strength. 112 and U4 are for- 
ward scattering interaction parameters. The boundary 
condition on the fermions at x = is determined by the 
transmission of the point contact, and will be discussed 
in various limits below. 



1. Bosonization of a single 



We first consider the Luttinger liquid theory for a sin- 
gle edge. We thus bosonize according to 

W = — L=e^, (2.3) 



\/27rx c 

where a = in, out, and x c is a short distance cutoff. ipi >a 
obey the Kac Moody commutation algebra, 

[<t>i,a{x),<t>3,b{v)] = inSijT^sgnix - y). (2.4) 
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into a theory of decoupled chiral bosons 
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where <p^ a obey 

[</>»,<, (x), 4>j,b(v)] = i7r 9 S v T ab s S n ( x - v)- 
Here v = vq^J (\ + A4) 2 — \\ and 



1 + A 4 - A 2 
1 + A 4 + A 2 



(2.7) 



(2.8) 



(2.9) 



The Luttinger liquid parameter g determines the power 
law exponents for various quantities. For instance the 
tunneling density of states scales as p{E) cx SW+VflO/ 2 - 1 . 



2. Mapping to Spinful Luttinger liquid 

Consider an open point contact in a Hall bar geometry 
with edge states on the top and bottom edges which con- 
tinuously connect leads 1 and 2 and leads 3 and 4. We 
then define left and right moving fields with spin f , J. as 

</>i?T = h.,ia{-x)6{-X) + (/)2,out{x)9(x) 

4>h\ = <p2; m {x)9{x) + <pi,out(-x)Q(-x) 

</>l T = h,in(x)e(x) + 4>i, out {-x)0(-x) (2.10) 

4>Rl = (j>i; m {-x)0{-x) + <j>3,out(x)6(x). 

It is then useful to define sum and difference fields as 
1 



(2.11) 



where a = R, L = +, — and a =T,I= +, — • Then, 9 a 
and (p a obey, 



[9 a (x),tpp(y)] = 2Tri8 al3 9(x - y), 
and J2J3 [231) become^ 

1 



(2.12) 



a—a.p 



g a (d x <p a ) 2 + —{d x a ) 2 
9a 



where 



9 P =9, 9a = 1/5, 



(2.13) 
(2.14) 
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TABLE I: The effect of discrete symmetry operations on the 
boson fieids 9 P and 6 a . 

and g and v are given in the previous section. 

It is useful to list the effect of symmetry operations 
on the charge-spin variables, because symmetries con- 
strain the allowed tunneling operators. Charge conser- 
vation leads to gauge invariance under the transforma- 
tion ip p — > tp p + 5 p . The conservation of spin S z leads 
to invariance under tp a — > f a +8 a . The effects of 
time reversal and mirror symmetries is shown in Table 
I. Time reversal symmetry is specified by the operation 
©■0a(70 -1 = io"4>a9- The mirror A4x interchanges leads 
14 <-> 23 while My interchanges leads (12 <-> 34). 

3. Four Terminal Conductance 

The central measurable quantity is the four terminal 
conductance, defined by 

/. Y^C.Aj- (2.15) 

3 

where Ii is the current flowing into lead i. Gy is in 
general characterize by 9 independent parameters. In 
Appendix A we introduce a convenient representation for 
these parameters, which simplifies the representation of 
symmetry constraints. Here we will summarize the key 
points of that analysis. 

The presence of both time reversal symmetry and spin 
conservation considerably simplifies the conductance. It 
is characterized by three independent conductances 

( Ix \ = ( Gxx Gxy\(v x \ , w) 

\Iy J V Gyx Gyy )\ Vy ) 

Here Ix = h + I a is the current flowing from left to right 
in Fig. 1, while Iy = I\ + I2 is the current flowing from 
top to bottom. Similarly, Vx is a voltage biasing leads 
(14) relative to (23) and Vy biases leads (12) relative to 
(34). Gxx is thus the two terminal conductance mea- 
sured horizontally, while Gyy is the two terminal con- 
ductance measured vertically. Gxy — Gyx is a "skew 
conductance" , which vanishes in the presence of mirror 
symmetry. Given these three parameters, the full four 
terminal conductance matrix Gjj can be constructed us- 
ing Eq. CM). 

A second consequence of spin conservation is the quan- 
tization of a particular combination of Gij. In particular, 
in appendix A we define a third current Iz = Ii + ^3 and 



a third voltage Vz which biases leads (13) relative to (24). 
Spin conservation then requires 

Iz = GzzVz, (2.17) 

with 

e 2 

G ZZ = 2-. (2.18) 

Since spin nonconserving spin orbit terms are allowed, 
spin conservation will not be generically present in the 
microscopic Hamiltonian of the junction. Nonetheless, 
we will argue that the low temperature fixed points pos- 
sess a emergent spin conservation, as well as mirror sym- 
metry, so that (|2.18p should hold, albeit with corrections 
which vanish as a power of temperature. 

B. Stable Phases 

In this section we describe various stable fixed points 
which admit simple descriptions using bosonization. We 
will first focus on the limit in which spin is conserved at 
the junction. There arc then four simple fixed point o 15 ' 16 . 
These include the perfectly transmitting (CC) limit, in 
which both charge and spin conduct, and the perfectly 
reflecting limit (II) in which both charge and spin are 
insulating. The "mixed" fixed points, denoted IC (CI) 
are perfectly reflecting for charge (spin) and perfectly 
transmitting for spin (charge). 

In the presence of spin non conserving spin orbit terms 
(which preserve time reversal symmetry) an additional 
fixed point is possible in which time reversal symmetry 
is spontaneously broken. We will see that in the presence 
of spin orbit terms this time reversal breaking insulator 
(TBI) phase is the stable phase when g < 1/2. 

1. Weak backscattering (CC) limit 

We first consider the limit where the point contact is 
nearly open and assume spin is conserved. It will prove 
useful to follow Ref. [TH and write l|2.13p as a 0+1 di- 
mensional Euclidean path integral for 6p_ a {r) = 9 p ^ a (x = 
0,t). This formulation is not essential for carrying out 
the perturbative analysis of this fixed point. However, 
it is of conceptual value for discussing the duality be- 
tween different phases, which can be understood in terms 
of instanton processes in which 9 Pj(7 (t) tunnels between 
degenerate minima at strong coupling. This is accom- 
plished by setting up the path integral for ap (x,T) and 
then integrating out 9 tT ^p(x,r) for x ^ 0. The resulting 
theory for CTi p(r) has the form of a quantum Brownian 
motion mode l ' 25 ' 28 ' 29 ' 30 , described by the Euclidean ac- 
tion 

1 1 r@ rj T 

S CC = TlY, KH*a("n)| 2 - / -Vcc(9 a ,9p), 

P £tf n 27r 9" Jo T c 

(2.19) 
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FIG. 2: Schematic representation of tunneling processes in 
(a,b,c) the CC phase ("small v") and (c,d,e) the II phase 
("small t"). (a) and (d) describe single electron processes, 
while the others are two particle processes. The duality re- 
t a and Vcr *-* t p can clearly be seen. 



where uj n — 2im/ (3 are Matsubara frequencies, and 
j3 = 1/ksT. We have included the short time cutoff 
t c = x c /v in the second term to make the potential 
V(9 pi 8 a ) dimensionless. The theory can be regularized 
by evaluating frequency sums with a exp(— |w„|r c ) con- 
vergence factor. 

The potential V(9 pi 9 (7 ) is given by an expansion in 
terms of tunneling operators, which represent the pro- 
cesses depicted in Fig. [2{a,b,c), 

Vac = v e cos(8 p + rjp) cos9 a + v p cos20 p + v a cos29„- 

(2.20) 

v e represents the elementary backscattering of a single 
electron across the point contact. The phase of cos# CT 
in that term is fixed by time reversal symmetry. The 
phase r\ p of cos 9 P is arbitrary, though mirror symme- 
try, if present, requires r\ p = mr. In addition we in- 
clude compound tunneling processes. v p represents the 
backscattering of a pair of electrons with opposite spins. 
We have chosen to define 6 p such that the phase of this 
term is zero. Note that this process involves the tun- 
neling of spin (not charge) between the top and bottom 
edges. Similarly, v a represents the transfer of a unit of 
spin from the right to the left moving channels, and in- 
volves the tunneling of charge 2e between the top and 
bottom edges. In general higher order terms could also 
be included. However, those terms are less relevant. 

The low energy stability of this fixed point is deter- 
mined by the scaling dimensions A(v a ) of the perturba- 
tions, which determine the leading order renormalization 
group flows, 



dv a /d£ = (1 — A(v a ))v a . 
These are given by 

A(«e) = (3p + £ ct )/2 = (.9 + .9" 1 )/2 



(2.21) 
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2g„ 



= 2.9 
= 2g 



It is therefore clear that all operators are irrelevant for 
1/2 < g < 2, so that the CC phase is stable. For g < 1/2 
v p becomes relevant, and for g > 2 v a becomes relevant. 



At the fixed point the conductance matrix elements are 

M t 

(2.23) 



Gxx = 2e 2 /h 
Gyy = Gxy = 0. 



At finite temperature, there will be corrections to these 
values. The leading corrections will depend on the least 
irrelevant operators. We find 



5G 



xx 



~ Cl v 2 e T3+9 l -a g > l/y/3 
9 < VV3 



-c 2 v 2 p T 4 9- 



6GYY ~ { c^T^- 2 g > V3 ' (2 ' 24) 

where Cj are nonuniversal constants. Note that for 
.9 < l/y/3 the exponents for Gxx and Gyy are dif- 
ferent. In addition, there will be power law corrections 
to Gxy when the mirror symmetries M. Xl M. y are vio- 
lated. However, this correction is zero when computed 
from (I2.19|2.20|) . even when r\ p ^ 0, due to the sym- 
metry of ()2.20j) under 9 a — ► — 9 a . Computing Gxy re- 
quires a higher order irrelevant operator. For instance 
\idxtptrSm9pCOS00 and \id x kp p cos9 p shi9 a - break both 
M.x and My, while preserving time reversal. This leads 
to 



SGxy = csXiXiT^ 



(2.25) 



Note that the temperature exponent of Gxy is at least 2 - 
even for weak interactions g ~ 1. This is because the tun- 
neling terms Ai and A2 include an extra derivative term. 
This is related to the fact (which we will show in Section 
III.C) that for non interacting electrons Gxy = 0. Weak 
interactions then introduce inelastic processes which give 
Gxy oc T 2 . The vanishing of Gxy is a unique property 
of the spin filtered edge states of the QSHI, which does 
not occur for a generic four terminal conductance. 



2. Weak Tunneling (II) limit 

When the point contact is pinched off, 6* Pi<T are effec- 
tively pinned, and a theory can be developed in terms of 
electron tunneling process across the point contact. This 
theory is most conveniently expressed in terms of the dis- 
continuity 9 a .p = (p" S p ht — <P&p across the junction 31 . The 
theory takes the form 



S 



11 



-y 



9a_ I 

2tt' 



9«(w»)| S 



dr 



Vu(o a ,e p ), 



'c 



(2.22) w ith 



(2.26) 



t p cos 29 p + t a cos 29 o 



Va = t e cos(9 p + rip) cos ( 

(2.27) 

As depicted in Fig. (^d,e,f) t e represents the tunneling 
of a single electron from left to right across the junction. 



7 



t a describes the transfer of a unit of spin across the junc- 
tion. t p describes the tunneling of a pair of electrons with 
opposite spins. 

The relationship between Sn and Sec can be under- 
stood in two ways. First, since both Sn and Sec de- 
scribe tunneling between the middles of two disconnected 
Luttinger liquids (either on the top and bottom of the 
junction or the left and right) the two theories are iden- 
tical. It is straightforward to see that if we make the 
identification 



2k 



(a) 



■A 



2k 



2k 



(b) 



2n 



(2.28) 



it follows that 



Sii{g P ,ga,t e ,tp,t,j) = Scc(g<j,gp,Ve,v,j,v p ). (2.29) 

Thus, the "small v" and "small t" theories are dual to 
each other, with the identification 

V e <-> t e 
Vn <-> t a 



Vn 



t 

<-+ g 



p 



(2.30) 



Using this identification, the scaling dimensions A(i Q ) 
can be read off from Eq. 12.221 Thus, like the CC phase, 
the II phase is stable when 1/2 < g < 2. The low tem- 
perature conductance can also be read from (|2.23[) , (|2.24|) 
and (|2.25|l using the identification 



Gxx Gyy- 



(2.31) 



Another way to understand this duality, which will 
prove useful below, is to consider an instanton expansion 
for strong coupling. For large v e (9 p ,9 a ) will be tightly 
bound at the minima of V(6 P ,6 <T ), shown in Fig. EJa). 
(Here we assume for simplicity r\ p = 0.) The partition 
function describing the path integral of (|2 . 19[) can then 
be expanded in instanton processes, in which (9 p , 9 a ) 
switches between nearby minima at discrete times. Eval- 
uating the first term in (|2.19|) for a configuration of in- 
stantons leads to an interaction between the instantons 
which depends logarithmically on time. The expansion 
describes the partition function for a one dimensional 
"Coulomb gas", where the "charges" correspond to the 
tunneling events. This Coulomb gas has exactly the same 
form as the expansion of (|2.26p in powers of t e , t p and 
t a . Thus, we can identify t e , t p and t a as the fugacity of 
the instantons. 

This duality argument also works in reverse. Start- 
ing from (|2.26[) we can derive (|2.19[) by considering large 
t e and expanding in instantons in 9 P and 9„ connecting 
minima in Fig. EJb), which have fugacities v e , v p and v a . 



3. Charge Insulator /Spin Conductor (IC) 

We next study the mixed charge insulator spin con- 
ductor phase. To generate the effective action for this 



49 r 



(c) 



2n 



A9r 



271 



-1 T z =1 



2k 



T x =-1 



(d) 



2k 



FIG. 3: (a) Positions of the minima of V(8 P , 9 a ) in Eq. [27201 
When the minima are deep instanton tunneling events be- 
tween the minima, denoted by t e , t p and t a correspond to 
the transfer of charge and spin across the junction, and de- 
fine the dual theory (|2.26I2.27[1 . (b) Positions of the minima 
of V(9 p ,Qcr) in the dual theory (12.26I2.27[) . Instanton process 
v e , v p and iv correspond to backscattering of charge and spin 
in the original theory, (c) The IC phase viewed from the CC 
limit. When v p is large and v a = 0, the minima of V(8 P , 9 a ) 
in Eq. 12.201 are on one dimensional valleys, and define the 
IC phase. When v„ is small but finite the valleys have a pe- 
riodic potential v a r z cos # CT , with opposite signs r 2 = ±1 in 
neighboring valleys. Instanton tunneling processes between 
the valleys, denoted t p , switch the sign of r z . (d) The IC 
phase viewed from the II limit, in which t p = and t a is 
large. The valleys have periodic potential t p r x cos 9 P with 
t x = ±1, whose sign is switched by instanton processes v a - 



phase, including the leading relevant operators it is use- 
ful to use the instanton analysis discussed at the end of 
the previous section. Consider (|2.26I2.27[) for large v p , 
keeping v e and v a small. 9 p will be pinned in the minima 
of — cos 26^, 8 p — n7r, while a remains free to fluctuate. 
{Opi^o) are thus confined to "valleys" along the vertical 
lines in Fig. [3](c) . 

There are two types of perturbations to be considered. 
First, v e will lead to a periodic potential along the verti- 
cal lines, with minima at the dots. Note, however, that on 
alternate lines the sign of the periodic potential changes, 
since cos 6 P cos8 a ~ (— l) n cos# CT for 9 p = nn. 

Next consider an instanton process where 8 p tunnels 
between neighboring valleys. In this process, 9 p — > 9 p ±n, 
but 9„ is unchanged. It follows that the v e perturbation 
discussed above changes sign. Thus, the instanton pro- 
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cess does not commute with the v e term. 

The expansion of the partition function in both instan- 
tons and v e can be generated by the action for the IC 
phase give by Sic = Sj C + Sj c with 



1 



ft 2tt 2ng, 



K||^K)| 2 , 

(2.32) 



and 



S 



ic 



dr 

r c 



t p T X COS Op 



V a T cos( 



(2.33) 



Here t p describes the instanton tunneling process. The 
tilde distinguishes it from the ordinary charge tunneling 
process, which involves charge 2e. t p describes a tunnel- 
ing of charge e without spin. v a describes the periodic 
potential as a function of 9 a generated by v e . We have 
introduced a pseudo spin degree of freedom t z = ±1 
to account for the sign of cos# p in the different valleys. 
Since the instanton process switches the sign, it is asso- 
ciated with t x . Expanding the partition function defined 
by (|2.32I2.33[) in powers of t p and v a precisely generates 
the expansion of (|2.19l2.20p in instantons. 

It is also instructive to derive (|2.32l2.33p starting from 
the opposite limit of the II phase described by (|2.26I2.27[) . 
In this case, consider large t a , which leads to the horizon- 
tal valleys as a function of 8 P and 9 a in Fig. G^d) . The 
roles of the two terms in (|2.33|) are thus reversed. t p de- 
scribes the periodic potential along the valleys, which has 
a sign specified by r x = ±1. v a describes the instanton 
processes which switch the sign of t x . 

The lowest order renormalization group flows depend 
only on the scaling dimensions of t p and v ai and are un- 
affected by the pseudospin t x,z . We find 



A(t P ) 
A(^) 



1 

9o 



1 

= % 
1 



(2.34) 



Thus, the IC phase is stable when g < 1/2. 

In section III.D we will require the renormalization 
group flow to third order in t p and v„. There, the non 
trivial interaction between them introduced by the pseu- 
dospin will play a crucial role. 

The conductivity at the IC fixed point is given by 



G 



xx 



G 



YY 



G 



XY 



0. 



(2.35) 



This, however, does not mean that the full four termi- 
nal conductance is zero because spin conservation still 
requires Gzz = 2e 2 /h. This leads to the non trivial 
structure in the four terminal conductance predicted in 
Ref. M 

At finite temperature, there will be corrections to the 
conductance. We find 



dit z p T3 



SGxx = 
SGyy = d 2 v 2 cr T^ 1 - 2 . 



(2.36) 



As in section the corrections to Gxy will depend 
on a higher order irrelevant operator. For instance, 
XiT y sin 9 p sin 9 a and \2T y cos 9 p cos 9 a lead to 



SGxy — d^Xi^T 29 



(2.37) 



As in (|2.25[) . Gxy is suppressed more strongly at low 
temperature than Gxx and Gyy, and the exponent is 
larger than 2 for g < 1/2. 



4- Charge conductor/Spin insulator (CI) 



For g > 2 the perturbation v a cos 29 a in (|2.20[) becomes 
relevant and drives the system to the CI phase. This may 
be described in a manner similar to the IC phase. It is 
described by the action Sci — S CI + S CI with 

1 



and 



S " = I°T 

JO T c 



2tt 



t„T Z C0S( 



UJn\\9p{0J n )Y 



(2.38) 



v p t x cos(9 p + rip) . 

(2.39) 



The leading relevant operators have dimensions 

2g a 2 



A(£ P ) = f = |- (2-40) 
This phase is thus stable when g > 2 and has conductance 



Gxx = G YY = 2e 2 /h 
Gxy = 0. 



(2.41) 



5. Spin orbit interactions, and the T-Breaking Insulator 

In this section we consider the role of spin orbit inter- 
action terms which violate the conservation of spin S z , 
but respect time reversal symmetry. We will argue that 
such terms are irrelevant for the critical behavior of the 
point contact when g > 1/2, but they are relevant for 
g < 1/2 and drive the system at low energy to a time 
reversal symmetry breaking insulator (TBI). 

Time reversal symmetry allows the following terms in 
the expansion about the CC fixed point (|2.19|) . 

Sec — / — [v so cos ip a smO a + v s f cos(2(p a + risf)] . 
Jo T c 

(2.42) 

The first term is a single electron process ^JjjV'-Rl (Fig. 
H£a)) in which an electron flips its spin and crosses the 
junction. The second term is a correlated tunneling pro- 
cess ijj^ip^ipniipLiiFig- HKb)), where a left and right 
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moving pair of up spins flip into a left and right mov- 
ing pair of down spins. Referring to Table HI it is clear 
that both terms respect time reversal symmetry. rj s f is 
allowed by time reversal symmetry, but violates both mir- 
rors M. x and Ai y . Higher order processes are also possi- 
ble, though they will be less relevant perturbatively. 

It is straightforward to determine the scaling dimen- 
sions of these perturbations. We find, 



A(u BO ) 



2 {9 ° 



1 



9a 



A(v sf ) = — =2. 9 . 

9a 



(2.43) 



For g ^ 1 the single particle spin orbit term, v so is always 
irrelevant. However, v s f becomes relevant when g < 1/2. 

At finite temperature these lead to corrections to the 
conductance of the CC phase. To lowest order they do 
not affect Gxx, Gxy and Gyy- However we find 



SGzz oc 




2 3>l/\/3 
9 < 1/V3 



(2.44) 



Like Gxy, Gzx and Gzy are zero unless higher order 
irrelevant operators, which involve extra powers of d x <p a 
or d x 9 a , are included. We find 



SGzx oc T 2 3 
SGzy oc T g+9 ~ 



(2.45) 



For weak interactions, g ~ 1 these conductances vanish 
forT^OasT 2 . 

For g < 1/2 there are two relevant perturbations about 
the CC limit. To study their effects we consider a model 
in which only the relevant perturbations appear. Since 
these perturbations involve the commuting operators <p a 
and 9p, it is useful to study the + 1 dimensional field 
theory of those variables 



E 



i 



2ir 9l 



KII^K)] 1 



2vr 



|w n ||<^ CT (w„)| 2 



with 



°cc 



(2.46) 



dr 



v p cos(26> p + rjp) + v s f cos(2ipo- + r) g f)] 



(2.47) 

The low temperature behavior of this theory can be stud- 
ied by the duality arguments of section II. B. 2. When v p 
and v s f are both large, (6 p , Lp a ) will be stuck in the deep 
minima of Vcc{@p, fa) shown in Fig. [5] In this phase, 
the four terminal conductance is zero, 

Gab = 0. (2.48) 
This can be seen most simply by renaming the variables 



"p 

fa 
fp 



6*i — 9 2 

fl + f2 
fl - f2- 



(2.49) 



(a) v so I or f 



(b) v sf 1 1 



FIG. 4: Tunneling processes in the CC limit allowed by spin 
nonconserving spin orbit interactions. v ao is a single parti- 
cle process where a single spin is flipped, while v s f is a two 
particle process, flipping two spins. 



271 



to 



9a 



A 



FIG. 5: 



Minima of the potential V{9 p . / tp (T ) in 12.471 Large 
v p and v so define the time reversal breaking insulating phase. 
Instanton processes t p and to- restore time reversal invariance. 
They correspond to tunneling of spinless chargons or charge- 
less spinons. 



The interpretation of #1(2) and fi(2) is simple. They are 
the usual Luttinger liquid charge and phase variables for 
the top (bottom) edges in Fig. [2la,b,c). In the strong 
coupling phase 9\ and #2 are both pinned, so that any 
current flowing in from any lead is perfectly reflected back 
into that lead. The four leads are completely decoupled. 

This is the same perfectly reflecting phase that would 
arise if we had a single particle backscattering term on 
each edge t>back(cos2#i + cos26> 2 ) = 2vbackCos0 p cos<^ tT , 
which would be relevant for g < 1. However in our prob- 
lem that term is forbidden by time reversal symmetry. 
It is thus clear that time reversal symmetry is violated 
by the strong coupling fixed point. It is useful to see 
this from Fig. [5] Note that since under time reversal 
(f a — > ip a + 7r. Thus pinning (p a violates time reversal. 
There are two sets of minima of V{9 p , ip a ) which are in- 
terchanged by the time reversal operation. 

At finite temperature tunneling processes between the 
two sets of minima of V(6 P , ipa-) will restore time reversal 
symmetry. These instanton processes correspond to tun- 
neling of charge from one lead to another. Interestingly, 
the lowest order instanton processes, denoted t p and t a 
do not correspond to tunneling of electrons, but rather 
spinless charge e "chargons" , or charge neutral "spinons" . 

The scaling dimensions of these instanton processes 
can be deduced from (|2.46I2.47|) . We find 



A(? P ) 



1 



1 

2^ 
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A(i.) = f = (2- 50 ) 

Thus, both processes are irrelevant for g < 1/2, and the 
TBI phase is stable. These processes lead to power law 
temperature behavior, 

SG XX = Cl P p T 1/s ~ 2 

5G Y Y = c^T 1 ' 9 - 2 . (2.51) 

When the t p ^„ processes dominate, there will be non 
trivial noise correlations in the current. The t p process 
involves transferring charge e/2 from lead 1 to lead 2 
and another e/2 from lead 4 to lead 3. This leads to 
correlations in the low frequency noise defined by 

SiiH= / ' dte iut {Ii{t)I j {0) + Ij(Q)Ii(t)). (2.52) 

Consider the two terminal geometry Ix — Gxx Vx ■ The 
current Ix will be carried by the t p processes, so that 
1 1 = I4 = Ix/2. The shot noise correlations in the limit 
lu — ► will be 

Sn = S44 = S*i4 = S41 = 2e* I\ (2.53) 

with e* = e/2. Thus, the currents are all perfectly corre- 
lated, and the current in each lead is carried by fractional 
charges, e/2. 

III. CRITICAL BEHAVIOR OF 
CONDUCTANCE 

In this section we describe the critical behavior of the 
conductance at the pinch-off transition of the point con- 
tact. We will compute the critical conductance G*, the 
critical exponent a g and the scaling function Q g {X) in 
certain solvable limits. We will begin in section IIIA 
with a discussion of the general properties of the scaling 
function and a summary of our calculated results. Then 
in the following sections we will describe in detail our 
calculations for g = 1 — e, g = l/\/3 and g = 1/2 + e. 

A. Scaling behavior and summary of results 

The stability analysis of the previous sections leads to 
the phase diagram as a function of g depicted in Fig. [S£a). 
The top line depicts the CC phase and the bottom line 
depicts the II phase, and the arrows denote the stability 
associated with the leading relevant operators. Since the 
II and CC phases are both stable for 1/2 < g < 2 they 
are separated by an intermediate unstable fixed point P, 
denoted by the dashed central line. For g < 1/2 the II 
and CC phases become unstable, and when spin is con- 
served the flow is toward the IC phase. We will see in sec- 
tion III.D that the unstable critical fixed point matches 
smoothly onto the IC fixed point at g — 1/2. Similarly, 



the CI fixed point is stable for g > 2, and connects to the 
critical fixed point at g = 2. 

For 1/2 < g < 2 the unstable intermediate fixed point 
P describes the critical behavior of the pinch-off transi- 
tion of the point contact. We will argue that this fixed 
point is characterized by a single relevant operator, which 
allows us to formulate a single parameter scaling theory 
for the pinch-off transition. If we denote u as the rel- 
evant operator, then the leading order renormalization 
group flow near the fixed point has the form, 

du/al£ — a g u, (3-1) 

where a g is a critical exponent to be determined. By 
varying a gate voltage Vq it is possible to cross from 
the region of stability of the II phase to the region of 
stability of the CC phase. In the process one must pass 
through the fixed point u = at Vg = Vq. Near the 
transition, we thus have u cx AVg — Vg — Vq. Under a 
renormalization group transformation in which energies 
length and time are rescaled by b, we have u — > ub a <> and 
T — > Tb. Invariance under this transformations requires 
that physical quantities can only depend on u and T in 
the combination u/T Q '. Close to the transition we thus 
have 

t,1^o Gab{T ' AVg) = 2 J^ab(c^), (3.2) 

where c is a nonuniversal constant and Q 9 ,ab is a univer- 
sal crossover scaling function which varies between and 
1. 

We will argue that the critical point characterizing the 
pinch-off transition has emergent spin conservation as 
well as mirror symmetry, so that the only nonzero el- 
ements of the conductance matrix are Gxx and Gyy. 
Moreover, the duality considerations discussed in section 
III.C require that G 3 ,yy{X) and Q g ,xx(X) are related, 
so that they are both determined by the same universal 
scaling function, 

G g ,xx(x) = g g (x), 

g g ,YY(x) = g g {-x). (3.3) 

The scaling function Q g {X) has some general proper- 
ties which are easy do deduce. First, the equivalence 
between the CC theory at g with the II theory at 1/g 
leads to the relation 

g 1/g (x) = i-g g (-x). (3.4) 

Second, when T — > for fixed AVg the system flows to 
either the CC or the II phase, where the temperature 
dependence of the conductance is given by (|2.24p . The 
behavior of the scaling function for large X then follows, 

g g {X^+oo) = l-a+X-ti 
g g {X^-oo) = a-X-K. (3.5) 
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FIG. 6: (a) Phase diagram for a point contact in a QSHI as 
a function of the Luttinger liquid parameter g. The arrows 
indicate the stability of the CC, II, CI and IC phases, as 
well as the critical fixed point P. This figure assumes spin 
conservation. In the presence of spin orbit interactions, the 
IC phase is unstable for g < 1/2. This leads to the TBI 
phase discussed in section II. B. 5. (b) Conductance G* of the 
critical fixed point as a function of g. The curve is a fit, 
which incorporates the data in l|3.7)l . (c) Critical exponent 



a g as a function of g. The curve is a fit incorporating the 



data in (|3.8[) . g is plotted on a log scale in all three panels to 
emphasize the g <-> 1/g symmetry. 



The coefficients a g depend on the normalization of X, 
but can be fixed if we specify G' g (X = 0) = 1/2. The 
exponents obey the relations 



# = 



_ J (4g-2)/a g 

(g + g- 1 -2)/a g 



= (g + g- 1 ~2)/c 



1/2 < g < 1/V3 
1/V3 < g < 1 

1/2 < g < 1. (3.6) 



The behavior of g for 1 < g < 2 can be deduced using 

§3). 



In the following sections we compute properties of the 
scaling function at g = 1 — e, g = l/v3 and g = 1/2 + e. 
From (|3.4[) we can deduce corresponding results at g = 
1 + e, g = V3 and g = 2 — e. First consider the critical 
G g (X = 0). We find, 



conductance G* 



1/2 + 0(e 3 ) .9 = l-e 
(V. 1 , : I (V3-l)/2 g=l/y/3 
7r 2 e g = l/2 + 



(3.7) 



These results are summarized in Fig. [6][b). The curve is 
a polynomial fit of G* (log g) which incorporates the data 
in Eq. (|3.7j) and the g <-> 1/g symmetry. It is satisfying 
that the curve is smooth and monotonic, which indicates 
a consistency between the slopes at g = 1/2, 1 and the 
value at g = l/y/3. 

We are able to deduce the critical exponent a g for g = 
1 - e and g = 1/2 + e. We find 



= I e 2 /2 5 = l-e 
1 4e g = 1/2 + 



(3-8) 



These results are summarized in Fig[6jc). The curve is 
a polynomial fit of a (log g). It is suggestive that in this 
fit a.g exhibits a maximum near g = l/v3 with a value 



h/VS 



.123 



1/8. It is possible, however, that a„ 



exhibits a cusp at g = l/\/3 analogous to the behavior 
of (5 g ii 



In sections III.C and III.D we compute the full scaling 



function Q g (X) in the limits g = 1 
lowest order in e. For g = 1 + e, e - 



e and g = 1/2 + e to 
* we find 



X 



^TTx 1 



o 



For g = 1/2 + e, e 

g 1/2 (X) = 0(X) 



1+X 



(3.9) 



(3.10) 



The singular behavior near X = in (|3.10[) is rounded 
for finite e. The perturbative analysis in Section III.D. 1 
shows that for \X\ <£L 1 



Gl/2+e(X) = 



X 



I _ e -X/(7rVT 



(3.11) 



Gi-t(X) and Gi/2+e{X) are plotted in Figs. El^a) and 
[7Jb). For g close to 1 the pinch-off curve is symmetri- 
cal about G* = e 2 /h. However, for stronger repulsive 
interactions it becomes asymmetrical, as G* is reduced, 
approaching at g = 1/2. 

The asymptotic \X\ — » oo behavior (|3.5[) of Gi(X) 
and Gi/2+e can also be determined from (|3.9l3.10p . 
though a separate calculation (see III.D. 3) is required 
for Gi/2+e(X — > — oo). The results, which are consistent 
with (|3.6[) are shown in Tabic HT1 
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A Gl-e(X) 




FIG. 7: The universal scaling function Q g {X) for (a) g — 1 — e 
(Eq. EU and (b) g = 1/2 + e (Eq. |330"]| . In (b) the solid line 
is e — > 0, and the dashed line shows the approximate behavior 
for e ~ .02. 



3 


# 


a+ 


&~ 


a g 


1 - e 


2 


1/4 


2 


1/4 


1/2 + e 


1 


1 


l/(8e) 


( 2 . 75) i/(8,) 



TABLE II: Parameters in Eq. I3.5l for the asymptotic behavior 
of the scaling function Q g (X) in the solvable limits g — > 1, 



B. Quantum Brownian Motion Model, Duality and 

s = i/Vs 

In this section we recast the Luttinger liquid model as 
a model of QBM in a periodic potential. This mapping 
elucidates the duality between the CC and II limits and 
exposes an extra symmetry the problem at g = 1/ y/3 
which allows us to deduce the critical conductance at 
that point. We begin with a brief review of the QBM 
model and then derive its consequences for the scaling 
function G g (X) and G* 



1. Quantum Brownian Motion Model 

The QBM model^^iSi was originally formulated as 
a theory of the motion of a heavy particle coupled to 
an Ohmic dissipative environment modeled as a set of 
Caldeira Leggett oscillators 3 ^. Though the applicability 



of this model to the motion of a real particle coupled to 
phonons or electron-hole pairs has been questione d 33 ' 34 , 
it was later shown that this model is directly relevant 
to quantum impurity problems. Specifically, the QBM 
model in a one dimensional periodic potential is equiv- 
alent to the theory of a weak link in a single channel 
Luttinger liqui d 16 ' 18 . In this mapping the QBM takes 
place in an abstract space where the "coordinate" of the 
"particle" is the number of electrons that have tunneled 
past the weak link. The periodic potential is due to the 
discreteness of the electron's charge. The low energy ex- 
citations of the Luttinger liquid play the role of the dissi- 
pative bath, and the strength of the dissipation is related 
to the Luttinger liquid parameter g. The one dimensional 
QBM model has two phases: a localized phase with con- 
ductance G = stable for g < 1 and a fully coherent 
phase with perfect conductance stable for g > 1. 

The SLL model corresponds to a QBM model in a two 
dimensional periodic potential, where the "coordinates" 
are the spin and charge variables 6> Pi(T . This model is 
richer than its one dimensional counterpart because it 
admits additional fixed points which are intermediate be- 
tween localized and perfect. These fixed points were first 
found in the Luttinger liquid mode]™,* 1 - 6 ", and later for- 
mulated in terms of the QBM2£. For certain values of g p 
and g a these intermediate fixed points are related to the 3 
channel Kondo problem^ and the 3 state Potts models^ 4 -. 
However, those limits are not directly applicable to the 
QSHI model, where g p = X/g a =9- We will show that 
when <7=1/ V3 the critical fixed point of the QSHI point 
contact corresponds to the intermediate point discussed 
in Ref. for a QBM model on a triangular lattice. 

To formulate the QBM model we begin with the action 
(|2.19p and define new rescaled variables, 

d a = ii\[2g a \r a . (3-12) 
Then (|2~pg|) takes the form, 

n v.j 

(3.13) 

The periodic potential is characterized by reciprocal lat- 
tice vectors G = mibi +ra2b2. The primitive reciprocal 
lattice vectors bi^ correspond to the single electron back 
scattering processes, and are given by 



27riG-r(r) 



bi 



1 



V2 



i 



V2 



(V9~ P ,-V^)- (3-14) 



The Fourier components of the periodic potential are 
Vhi = v h2 = v e e" ) "/4, w bl +b 2 = v p /2 and v hl - h2 = v a /2. 

The dual theory is obtained by expanding the partition 
function for large vg in powers of instantons. When vg 
is large, the potential has minima on a real space lattice 
R = niai + U2B.2- The primitive lattice vectors satisfy 



ij and are given by 
, 1 1 , 



ai 



V2 ^/Tp 1 y/!fc' 



a 2 



\/2 V5p' V9°- 



). (3.15) 
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FIG. 8: (a) Minimia of the periodic potential V(r) in (I3,13|l . 
(b) Minima of V(k) in the dual theory (pTTo) . When g a = 3g p 
both periodic potentials have triangular symmetry at the crit- 
ical point, which implies the mobility /x*^ is isotropic. This 
occurs at g = l/\/3. 

The expansion in instantons connecting these minima is 
generated by the action 

^^Ei-«iikK)i 2 -E/f^ 2 " R ' k(r) - 

(3-16) 

This is equivalent to ^2.2612.27)) with k a = ny/g a /29 a 
and £ ai — t a2 — t e e ^ p /4, £ ai + a2 t p /2, £ ai _ a2 = t (T /2. 

With the above normalizations for r and k the scaling 
dimensions of the potential perturbations are 

A(v G ) = \G\ 2 ; A(t R )HR| 2 . (3-17) 

Since operators are relevant when A < 1, the most rel- 
evant potentials are those with the smallest lattice (re- 
ciprocal lattice) vectors |R m in| (|G m i n |). As shown in 
Refs. [l6l and [25l there are ranges of g p and g a where both 
|R m i n | and |G min | > 1, so that both phases are perturba- 
tively stable. An unstable intermediate fixed point must 
therefore be present between them. 

This fixed point can be accessed perturbatively when 
|R m in| an d |G min | are close to 1. While this does not oc- 
cur in the regime g p = l/g a relevant to the QSHI prob- 
lem, it is instructive to study this perturbation theory 



because it provides evidence that the critical fixed point 
has emergent mirror and spin conservation symmetry. 

When g p = 1/2 + e p and g a = 3/2 + e CT the period 
potential has triangular symmetry, which is slightly dis- 
torted if e a ^ 3e p . If we denote the relevant variables as 
vi = Vh 1 = Vh 2 = v e e lrip /4: and V2 = Vb 1 +b 2 = v p /2, the 
second order renormalization group flow equations are 16 

dvijdl = -(e p + e a )vi - 2v^V2 

dv 2 /d£ = 2e p v 2 - 2v\. (3.18) 

These equations describe an intermediate fixed point 
with a single unstable direction at v\ = \] t P {t P + e<r)/2 
and V2 = (e p + e<r)/4. Note that at the critical point v\ 
is real, so that r\ p = 0. Thus the critical point has an 
emergent mirror symmetry even if the bare parameters 
in the model do not. Moreover, the flow out of the fixed 
point along the single unstable direction is also along a 
line with v\ real. Thus the crossover between the in- 
termediate fixed point and the trivial fixed point, which 
determines the crossover scaling function also has emer- 
gent mirror symmetry. Mirror symmetry breaking is an 
irrelevant perturbation at the critical fixed point. 

If e CT = 3e p then the lattice vectors have a triangular 
symmetry. In this case, the fixed point is at v± — v 2 = e p - 
This means that the periodic potential at the fixed point 
has emergent triangular symmetry, even when the bare 
potential does not. The unstable flow out of the fixed 
point is also along the high symmetry line v\ = v 2 - 

It seems quite likely that the critical fixed point and 
unstable flows connecting it to the trivial fixed points 
retain their high symmetry even outside the perturbative 
small e regime. This suggests that in general the critical 
fixed point has mirror symmetry, and that at g = l/v3 
it has triangular symmetry. We will use this fact below 
to determine the critical conductance at g = 1/ a/3. 

2. Kubo conductance, mobility and duality relations 

The spin and charge conductances in the Luttinger liq- 
uid model computed by the Kubo formula are given by 
a retarded current-current correlation function. For the 
present discussion it is useful to write this as an imagi- 
nary time correlation function, which can be analytically 
continued to real time via iu> — > lo + in before taking the 
lu —> limit. Then 

G%(M = -^nj dre^{J a (r)J (O)}, (3.19) 

where the spin and charge currents are J a = edtOa/i: — 
e[6 ,H] I (iirh) . This may be expressed as 

e 2 

Ga^n) = 2— v /g a g f3 fi a/3 (uj n ), (3.20) 
where the mobility of the QBM model is 

Ma/3(^n) = 2Tr\uj n \(r a (-uj n )rp(uj n )). (3.21) 
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/i Q /3 is normalized so that when vq, = \L a p = 8 a p. 

The conductance - or equivalently [i a p can also be 
computed from the dual model. It is given by 



Ma/3 — <W ^ P>a0> 

where the dual mobility is 

AV>(w„) = 2ir\u n \(k a (-u n )kf3(uj n )) . 



(3.22) 



(3.23) 



(|3.22I3.23|) are obvious in the perfectly transmitting and 
perfectly reflecting limits. They can be derived more 
generally by starting with a Hamiltonian formulation of 
the action, analogous to (|2.13j) , which involves both r and 
k. n a /3 can then be computed either by first integrating 
out k to obtain p.21[) or first integrating out r to obtain 
(I3~23l) . 

Since g p = l/g a = g, the dual theory depicted in Fig. 
[5](b) is identical to the original theory shown in Fig. [SJa) 



with the identification r p «-> k„ 



the mobility n* a g of the fixed point satisfies 



k p . It follows that 
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(3.24) 



In addition, if u parameterizes the relevant direction at 
the critical fixed point, then under the duality u — > — u. 
It follows that slightly away from the critical fixed point 
we have 



(3.25) 



Properties (|3.22p and (|3.25p imply that fi pp (u) = 1 — 
M<t<t (-«)• Using (|3.2l3.20IA15p . this leads directly to the 
property (|3.3|) of the crossover scaling function. 

An additional set of relations follows from the equiva- 
lence between the theory characterized by g and the dual 
theory characterized by 1/g. From this we conclude that 

M S ,a/3( u ) = Al/ff.o^O")- (3-26) 

This, combined with lj3.2|3.2QI3.22IA15)l ) leads to ([3^4]) . 



It is well known that the physical conductance 
measured with leads is not given by the Kubo 
conductanc o 35 i 36 ' 37 i 38 ' 39 . Rather, the Kubo conductance 
needs to be modified to account for the contact resistance 
between the Luttinger liquid and the leads. In appendix 
A we review the relation between the physical four termi- 
nal conductance and the Kubo conductance. From (|A19|) 
we conclude that 



Gxx = G YY = (\/3-l)-. 



C. Weak interactions : g — 1 — e 



(3.30) 



In this section we develop a perturbative expansion 
for weak interactions to compute exactly the crossover 
scaling function Q g (X) as well as the critical exponent 
a g for g = 1 — e. A similar approach was employed 
by Matveev, Yue and Glazmar*2£ to compute the scal- 
ing function for the crossover between the weak barrier 
and strong barrier limits in a single channel Luttinger 
liquid. In the single channel problem the transmission 
for non interacting electrons is characterized by a trans- 
mission probability T. Weak forward scattering interac- 
tions lead to an exchange correction to T at first order 
in the interactions. This correction diverges for E — > Ep 
as log IE 1 — Ef\. Matveev, Yue and Glazmar*2£ used a 
renormalization group argument to sum the log divergent 
corrections to all orders, to obtain the exact transmission 
T(E). 

For non interacting electrons, the QSHI point contact 
is characterized by a 4 x 4 scattering matrix which 
relates the incoming wave in lead i to the outgoing wave 
in lead j, 



(3.31) 



In terms of Sa the four terminal conductance is 



3. Conductance at g = l/v3. 

When g = 1/ \/3 the lattice generated by bi and b2 has 
triangular symmetry. In section III.B.l we argued that 
this means that at the critical fixed point the periodic 
potential also has triangular symmetry. The Cg rota- 
tional symmetry of the triangular lattice requires that 
the mobility is isotropic: 

Ha/3 = MO<5q/3- (3.27) 

Combining (pT22]) . ([3~24| . and (|3~27| requires that 



1 



(3.28) 



It follows from (|3.20p that the Kubo formula spin and 
charge conductances are given by 

«S = V3i (3.29) 



Gij — , (Sij 



(3.32) 



Under time reversals | ^ iOUt(in )) = +(-)Q«|'0j,m(out)), 
where Q = diag(l, — 1, 1, — 1). This leads to the con- 
straint S — ~QS T Q. This combined with unitarity 
S^S = 1 allows S to be parameterized as 



( o 



S = £/ f 



V 



t 

t 

-/ r* 
r f* 



f r \ 

r * —f* 

-t* 

-t* / 



u. 



(3.33) 



where Uij = Sije lXi is an unimportant gauge transfor- 
mation. The complex numbers t and r describe the am- 
plitudes for spin conserving transmission and reflection 
across the point contact, while / describes the ampli- 
tude for tunneling across the junction, combined with a 
spin flip. / = if spin is conserved. The conductance 
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can be expressed in terms of the transmission probabil- 
ities TZ = \r\ 2 , T = \t\ 2 and T = |/| 2 , which satisfy 
TZ + T + T =\. We find 



Gxx 


2e 2 




Gyy 


2e 2 


(3.34) 


Gzz 






Gab 


= for A ^ B. 





For a generic four terminal conductance device time 
reversal symmetry guarantees only the reciprocity 
relation 42 Gy = Gji, (or equivalently Gab = Gba)- For 
the QSHI point contact, the spin filtered nature of the 
edge states leads to additional constraints. First, the am- 
plitude for an electron to be reflected back into the lead 
it came from is Sa = 0. Thus Gn = e 2 /h. A second 
less obvious constraint is that G13 = G24, which when 
combined with reciprocity and unitarity is equivalent to 
G12 = G34 and G14 = G23. This leads to the vanishing 
of the skew conductance Gxy as well as Gxz and Gyz 
even when mirror symmetries Aix and M.y are explic- 
itly violated. This is a property of the non interacting 
electron model and can be violated with electron electron 
interactions if the mirror symmetries are absent. 

In order to compute the renormalization of the S ma- 
trix due to interactions it is useful to study the pertur- 
bative expansion of the single electron thermal Green's 
function, which can be represented as a matrix in the 
lead indices i, j as well as the channel labels a = in/out. 



-1 ah 



(x,t;x',t') = -i{T T [^ a {x,T)i> Ux',t')]), (3.35) 



where T T denotes imaginary time ordering. For non in- 
teracting electrons we have 



<&ij (z, z') 



1 

27ri 



Sij 


rr* 
ji 


z—z' 


z—z' 


Sij 


Sij 


z—z' 


z—z' 



(3.36) 



where z — r + ix and z — t — ix, and the a =in/out 
indices are displayed in matrix form. 

We now compute the perturbative corrections to 
G° ut,ln using the standard diagrammatic technique. For 
simplicity, we adopt a model in which 114 — 0, so that 
the only interaction term involves w 2 (V'/ n V'in)(V ; outV'out)- 
This considerably simplifies the analysis because many 
of the diagrams are zero. For instance, the exchange 
diagram shown in Fig. [9ja), which was responsible for 
the renormalization in the single channel Luttinger liq- 
uid problem is zero because it must involve G)!^ 011 *. This 
off diagonal Green's function depends on Skk which is 
zero due to the time reversal symmetry constraint. From 
(EH), g = y/(2wv F - \ 2 )/(2irv F + \ 2 ) ~ 1 - X 2 /(2ttv f ). 
Thus for g = 1 — e we may replace u 2 by 2ttv f ^- The 
nonzero diagrams at second order in u 2 are shown in 



(a) 



j.out 



kjn k.out 



(b) 
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FIG. 9: Feynman diagrams for the single electron 
Green's function. The dashed line is the interaction 
M2(Vm^/'in)(V , lutV , out). The exchange diagram (a) vanishes be- 
cause it involves Skk, and diagrams (b) and (c) cancel one 
another, (d) and (e) lead to a logarithmic correction to the S 
matrix given in (|3.38p . 



Fig. IHKb-e). Evaluating the second order diagrams gives 
a Green's function of the form 



1 



2tti z — z' 



(3.37) 



with 



SijSjiS^ 



SikSkiSi k SkiSij 

kl 



(3.38) 

where A are E are ultraviolet and infrared cutoffs re- 
spectively. The first term in the brackets was due to the 
diagram in Fig. [DJd), while the second term was from 
Fig. (He). Diagrams Mb) and(5](c) cancelled each other. 
Rescaling the cutoff A — * Ae~ £ leads to a renormalization 
group flow equation for Sij, 



dSij 
~~d£ 



kl 



(3.39) 



It is useful to rewrite this in terms of the transmission 
probabilities T, TZ, T . The renormalization group flow 
equation then can be written in the form 



dT/dl = e 2 T{T -T 2 -TZ 2 - T 2 ) 
dK/d£ = e 2 K{K - T 2 - K 2 - T 2 ) 
dTjdl = e 2 T{T -T 2 -Tl 2 - T 2 ). 



(3.40) 



The flow diagram as a function of TZ, T and T is shown 
in Fig. [TUl There are seven fixed points. The bottom 
corners of the triangle are the stable fixed points at TZ = 
1, T = T = (the II phase) and T = 1, TZ = T = (the 
CC phase). The third stable fixed point at the top of the 
triangle with T = 1, T = TZ = 0, corresponds to the case 
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F=1 




R=1 



T=R=1/2 



T=1 



FIG. 10: Renormalization group flow diagram for the trans- 
mission probabilities T, 1Z and T based on (|3.40[) represented 
in a ternary plot. The CC, II and P fixed points of interest 
in this paper, which have T — are on the bottom of the 
triangle. 



where an incident electron is transmitted perfectly with a 
spin flip. This is presumably difficult to access physically. 
On the edges of the triangle are unstable fixed points 
describing transitions between the different stable phases. 
The critical fixed point P of interest in this paper is the 
one on the bottom of the triangle at7£ = T= l/2,.F = 
0. Note that at this fixed point the spin non conserving 
spin orbit processes, represented by T are irrelevant. At 
the center of the triangle, at 1Z = T = T = 1/3 is an 
unstable fixed point describing a multicritical point. 

To describe the critical fixed point P and the crossover 
to the II and CC phases we now specialize to T = 
and consider the flow equation for the single parameter 
T characterizing the point contact, 



dT/dl = -e 2 T(l - T)(l - 2T). 



(3.41) 



Eq. 13.411 can be integrated to determine the crossover 
scaling function. If at I = T = T°, then, 



i 



T° - 1/2 



-e 2 i 



y/(T°- 1/2) 2 +T (l-r°)e- 

(3.42) 

As the gate voltage Vg is adjusted through the pinch- 
off transition, T° passes through 1/2 at Vg = V G , so 
T° — 1/2 cx AVg- At temperature T we cut off the renor- 
malization group flow at Ke~ l cx T. The conductance 
is then given by G X x = 2(e 2 /h)T(£ = log(A/T)). For 
AV G , T^Owe define X = (2T° - l)e^ 1 / 2 cx AV G jT^I 2 
and write the conductance in the scaling form, 



G X x(AV G ,T) = 2—g 1 - 
h 



(3.43) 



T a 3 

where c is a non universal constant, the critical exponent 
is 



«!_ £ = e 2 /2, 



and 



X 



Vl + A 2 



(3.44) 
(3.45) 



We find that the logarithmic renormalization to the S 
matrix accounts for the only correction to the conduc- 
tance to linear order in e. In principle one must consider 
a "RPA like" diagram for the conductance evaluated by 
the Kubo formula. While this gives a correction for an 
infinite Luttinger liquid at finite frequency, the correction 
is zero for a finite Luttinger liquid connected to leads in 
the w — ► limi t 35 ' 36 i 37 ' 38 . Since the critical conductance 
satisfies G* = 1-G*, it follows that G\_ t = l/2+0(e 3 ). 



D. 



3 = 1/2 + , 



g = 1/2 is at the boundary where the CC and II 
phases become unstable and the IC phase becomes sta- 
ble. We will show that when g = 1/2 + e the critical 
fixed point describing the transition between the CC and 
II phases approaches the IC fixed point and can be ac- 
cessed perturbatively using theory developed in Section 
II. B. 3. In addition, when g = 1/2, the marginal opera- 
tors v p cos 26 p at the CC fixed point and t cos 9 p at the IC 
fixed point can be expressed in terms of fictitious fermion 
operators. This fermionization process allows the entire 
crossover between the CC and IC phases to be described 
using a non interacting fermion Hamiltonian. A simi- 
lar fermionization procedure can be used to describe the 
crossover between the II and IC phases, which connect 
the marginal operators v a cos Be and t a cos 26 a . This will 
allow us to compute the full crossover scaling function 
G g (X) for g = 1/2 + e. 

We will begin by discussing the perturbative analy- 
sis of the IC fixed point and then go on to describe the 
fermionization procedure. 



Perturbative Analysis 



The IC fixed point is described by (j232l l2"U3"j) . When 
g = 1/2 + e the perturbations t p r x cos9 p and t a T z cosd a 
both have scaling dimension A = 1 — 2e, so the IC fixed 
point is weakly unstable. When v a = 0, nonzero t p is 
expected to drive the system to the CC phase, while for 
t p = nonzero v a will drive the system to the II phase. 
Thus, when both t p and v„ are non zero there must be 
an unstable fixed point which separates the two alterna- 
tives. This fixed point can be described by considering 
the renormalization group flow equations to third order 
in v a and t p . 

The first order renormalization group equation for t p 
is determined by the scaling dimension A(t p ). The next 
nonzero term occurs at order t p v 2 . To compute this term 
it is sufficient to use the theory at e = 0. Consider the 
third order term in the cumulant expansion of the par- 
tition function, when fast degrees of freedom integrated 
out: 

\ J dr 1 dT 2 {(T r [O p (r)O cr (T 1 )0 CT (T 2 )]) > 
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-<O p (r)) > (T T [0 <r (ri)O ff (7 S )]) > }. (3.46) 

Here O p = {I p /t c )t x cos 9 p and O a = {v a /t c )t z cos 8 a . 
T T indicates time ordering, and (•)> denotes a trace over 
degrees of freedom with A/6 < u> < A, and we assume 
for simplicity 6 > 1. Since # p and are independent 
and commute with one another the other disconnected 
terms all cancel. Moreover, the two terms in (|3.46|) will 
cancel each other unless the time ordering of the t x and 
t z operators leads to a relative minus sign between them, 

{T T [O p (T)O a (ri)O a (T2)})> = 

s±{0 P {T)) > {T T [O a {T 1 )O a (T 2 )]) > , (3.47) 

where s± = sgn(r — ti)(t — r 2 ). Thus the pseudospin 
operators in (|2.33[) play a crucial role in the renormal- 
ization of t p . Using the fact that (T t [0 (T (ti)0 (X (t2)])> = 
u^/2(ri — T2) 2 for e = we find that the third order 
correction to t p is 5t p = —t p v 2 \ogb. This leads to the 
renormalization group flow equation for t p , along with a 
corresponding equation for v a , 



dtp/dt = 2et p -t p vl 
dv a jdt — 2ev a — v a P 



(3.48) 



The renormalization group flow diagram is shown in Fig. 
1111 There is an unstable fixed point P at t p = v a = s/2e, 
with a single relevant operator. P separates the flows 
to the CC and II phases for which which t p or v a grow. 
Note that spin orbit terms such as v so and v s f discussed 
in Section II. B. 5 are irrelevant at P (see Eq. I2.43|) . This 
perturbative calculation provides further evidence that P 
exhibits emergent spin conservation, as well as emergent 
mirror symmetry. The critical exponent associate with 
the single relevant operator a P is 



«l/2+, 



= 4e. 



(3.49) 



The Kubo conductance G pp at the fixed point can be 
computed from (|3 . 1 9[) by identifying the current operator 



(t P /Tc) sin P . 



(3.50) 




FIG. 11: Renormalization group flow diagram characterizing 
the critical fixed point P for g — 1/2 + e. When v a and t p 
are small, the flows are given by (|3.48|l . On the axis v a = 
the fermionization procedure outlined in section III.D.2 de- 
termines the entire crossover between the IC and CC fixed 
points. A similar theory describes the crossover between the 



IC and II fixed points for t p 



0. 



where TZ a — tt 2 v 2 can similarly be interpreted as a re- 
flection probability for a different fictitious free fermion 
at g = 1/2. At the critical fixed point T p = lZ a = 2ir 2 e. 
Thus, 



G* 
XX 

G* 
XY 



G* 
YY 

0. 



y 2 

2— 7T 2 e. 
h 



(3.54) 



The behavior away from the critical point can be de- 
termined by integrating (|3.48[) . To this end it is helpful 
to rewrite (|3.48| in terms of T p and lZ a in the form 



d{T p - 71a)/ dl 
dlog(T p /TZ a )/d£ 

If (T p ,K a ) = (T^TZl) for 



4e(T p - K a ) 
(2/n 2 )(T p -U a ). 

= 0, then we find 



(3.55) 



This leads to 



(3.51) 



It is useful to define T p = n 2 t 2 p . We will see in the follow- 
ing section that this can be interpreted as a transmission 
probability for fictitious free fermions that describe the 
problem at g = 1/2. In terms of T p (noting that T p <C 1 
in this perturbative regime) we may use (|Al9j) to write 
the physical conductance as 



Gxx 

A similar calculation gives 



G 



YY 



h 



(3.52) 



(3.53) 



Ha(£) 



(T p °-K)e 



id 1 



exp 



0^„4e 



1 



fio 6Xp 



27r 2 e 



(e 4 ^ - 1) 



(3.56) 



At the pinch-off transition Vq = Vq, TZq = T . Thus, 
7q — IZo cx AVq- At temperature T we cut off the renor- 
malization group flow at Ae~ e oc T. Thus, in the limit 
AV G , T-^0we define X = {T p °-K° a )e 4e£ /2 oc AV G /T 4€ . 
The conductance then has the form 



G XX (AV G ,T) = 2-G g 
Gyy(AV g ,T) = 2jG g 



AV G 
c AV G 



(3.57) 
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with 



Gl/2+e(X) = 



X 



(3.58) 



This perturbative calculation is only valid when 
T p ,lZ a <C 1. Thus (|3.58p breaks down at low temper- 
ature, since as the energy is lowered either T p or 
grows. (|3.58|) is valid as long as \X\ <C 1. Note, how- 
ever that when e <C 1 we have Gi/2+e(X) = X9(X) when 
e <C X -c 1. In this regime, the smaller of T and 1Z has 
gone to zero. Thus we have 



4c f 



K a {£) = 
T p (£) = 



(T°-7e°)>0 



(K - T P °) > 0, (3.59) 



and the unstable flow is either on the x or y axis of Fig. 
rm In the next section we will solve the crossover ex- 
actly on these lines. This will allow us to compute the 
Gi/2+e{X) exactly for all X. 



2. Fermionization 

In this subsection we study the crossover between 
the IC fixed point and the CC and II fixed points for 
g = 1/2 + e. There are two cases to consider. First, for 
AVg > we will study the crossover between the IC and 
CC on the horizontal axis of Fig. [II] with v a = 0. This 
problem can be mapped to a single channel one dimen- 
sional fermi gas with weak electron electron interactions 
proportional to e. This allows us to use the method of 
Matveev, Yue and Glazmar*2£ to compute the crossover 
scaling functions G g ,xx{X) and G g ,YY{X) for X > ex- 
actly. For AVg < the crossover between the IC and II 
fixed points is on the vertical axis of Fig. [TTlwith t p = 0. 
This can be fermionized by introducing a different set of 
free fermions to compute the scaling functions for X < 0. 
The latter calculation (which is virtually identical to the 
former) is unnecessary, however, because we can use (|3.3|) 
to deduce the scaling functions for X < 0. We will there- 
fore focus on the IC to CC crossover. 

The crossover between the IC and the CC fixed points 
can be described by the action in the CC limit 



S cc=o>,„ — \u n \\0 p (u n )\ 2 + / dT Vp cos29 P 
P „ 2lr 9 J 



(3.60) 



v p <C 1 describes the CC phase. When v p ^> 1 the dual 
theory, formulated as in section II. B. 3 in terms of instan- 
tons with amplitude t p , describes the IC phase. When 
v a = at the IC fixed point we can safely ignore the 
pseudospin, and set t x = 1. 

For (7=1/2 this model is equivalent to the bosonized 
representation of a weak link in a single channel non in- 
teracting fermion with weak backscattering. 



Hf = -ivtjj^d x <7 z -ijj + Vfffl <j x i>5(x) 



(3.61) 



where ip = (ipu, iPl) T is a two component fermion opera- 
tor describing right and left movers. Using the bosoniza- 
tion relation (|2.3[) we identify 29 p — <j>R — <j>L and t>/ = 
nVp/v. The free fermion problem is solvable and charac- 
terized by a transmission probability T p = sech 2 (v / / v) . 
The free fermion solution therefore connects the CC limit 
(T p = 1) with the IC limit (% = 0). 

The Kubo conductance G p may be computed with the 

identification J p — dtOp/ix — v^a z ^p, giving 



G 



K 



h v 



(3.62) 



Note that this is the same as (|3.52j) . derived in the op- 
posite limit near the IC fixed point. When v p is large, 
T p C 1, and we can identify T p = (irt p ) 2 . The physi- 
cal conductance, measured with leads can be determined 
following the analysis in appendix A. From (|A19|) we find 



Gxx — 2- 



h2-T n 



Since v a = in (|3.60p . we have 



G 



YY 



= 0. 



(3.63) 



(3.64) 



For g = 1/2 + e the IC fixed point becomes slightly 
unstable, while the CC fixed point becomes slightly sta- 
ble. In this case the free fermion problem includes a weak 
attractive interaction 



i_/int 



(3.65) 



with Uf — 2ixve. This leads to a logarithmic renormaliza- 
tion of T p , which drives a crossover to the CC limit. The 
correction to T p occurs at first order in uf, and is due to 
the exchange diagram, shown in Fig. OJa). The analysis 
is exactly the same as that performed by Matveev, Yue 
and Glazman. As in section III.C the result can be cast 
in terms of a renormalization group flow equation for T p . 



dTp/d£ 

Integrating (|3.66|) gives 
T P {1) 



4eT p (l - T p ). 



T°e 



c 4ef 



4e£ 



1) ' 



(3.66) 



(3.67) 



where T p ° = T p (£ = 0). The scaling function for AVg > 
then follows by using the initial condition from (|3.59p . 
so that T p ° cx AV" G . Then, for AVg, T -> we define 
X = Tp°e 4e 72 oc AVb/T 4£ . Using (I3~B5D . (I3~MD and 
(|3.67p . the conductance has the scaling form for X > 

X 



GxX,l/2+e{X) 
GYY,l/2+e(X) 



X 

0. 



1 



(3.68) 



Using (|3.3|) . we may deduce the corresponding behavior 
for AVg < (or X < 0). The scaling function then has 
the form 



G 1/2+£ (X) = 6(X) 



X 



X + l 



(3.69) 
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Note that for X < 1 G\/%+ € {X) = Xd(X), in agreement 
with the limiting behavior of (|3.58[) for X ^> e. These 
two expressions can thus be combined to give 



X 



X + l- e- x /^*) ' 



(3.70) 



which reproduces (|3 . 58[) when \X\ ~ e <C 1 and (|3.69[) 
when |X| ^> e. This function is plotted in Fig. [TJb). 
Note, however, that this formula does not correctly cap- 
ture the leading behavior for X < when |X| ^> e. 
In particular, it misses the X — > — oo behavior, which 
([53]) and dSU) predict is proportional to |j^|-i/(8e) _ Xms 
regime is analyzed in the following section. 



3. Rebosonization 

We now analyze the leading behavior of G\/2+ € {X) for 
X < and |X| 3> e when e is small. Equivalently, we 
consider Gi/2+e,YY{X) for X > 0. This requires extend- 
ing the renormalization group flow equation for v a given 
in (|3.48|) to all t p (or equivalently T p ). This can be done 
by using the fermionized representation of t p T x cos 9 P in 
(I2.33p . The key point is that the presence of the pseu- 
dospin operator t x means that the operator v a T z cos 
changes the sign of the transmission amplitude for the 
fcrmions ip. This results in an X ray edge like contri- 
bution to the renormalization of v a . This can be com- 
puted by a method analogous to that used by Schottc 
and Schotte^i to solve the X ray edge problem, which in- 
volves transforming the non interacting fermions to even 
and odd parity scattering states and then rebosonizing. 
This approach was used to study the X ray edge problem 
in a Luttinger liquid in Ref. |41| . 

We begin by writing f2T33")) . H = H a +H p with 



nj -t_/0 



VrrT COS( 



and 



H p = -iv^a z d x rp + t f T x ^a x ^5{x). 



(3.72) 



Here Tt ° is the a part of (|2.13p , and we explicitly account 
for the pseudospin t x . Eq. 13.721 can be rebosonized by 
first replacing ip2(x) —>■ tp2(—x), which transforms the 
non chiral fermions to chiral fcrmions, eliminating the 
<j z in the first term, but leaving the second term alone. 
Then we perform a SU(2) rotation (^1,^2) — ► (V'ejV'o)) 
which changes a x in the second term into a z . ip e i ) de- 
scribe the even (odd) parity scattering states character- 
ized by scattering phase shifts 5 e = —8 that specify 
^Pe(o)( x > 0) = e 2 '" 5 "' ' ip e ( ){x < 0). We next bosonize 



"He 



V 

8tt 



•° /^/2ttx c and define <f>± = <f) e ± <f> . Then 



2n 



5^T x {d x( j>.)5{x), 



where 4>± obey, [4>±(x),cj>±(x')] = 27risgn(x — x'). 5. 
5 e — S a is related to the transmission probability by 



T p = sin 5- . 



(3.74) 



(5_ can be eliminated from (|3.73|) by the canonical 
transformation U — exp[iT x 5-<fi-(x — 0)/(27r)], which 
shifts 0_ —> (j>- + sgn(x)S-T x . This transformation also 
rotates r z in (|3.7ip , which becomes 



TC 



T+e*+- s -/*+T e 

(3-75) 

where r = t z ± ir y . The renormalization of v a can then 
easily be determined for arbitrary 5— . We find 



(3.76) 



For small t p , 6- = nt p , and (|3.T6[) reproduces (|3.48l) . 




However, (|3.76p remains valid to lowest order in e for all 

We now integrate (|3.55|) to a scale £q where from (|3.56p 
T p (e ) = 2X and K a {Zo) = 2X G e~ Xo '^ 2 ^ is small. 
(Here X = (T p ° - TZ a )e 4ee °/2.) We then use that as 
an initial value for (|3.76p , which we integrate assuming 
T p (£) is given by (|3.67|) and is unaffected by the small 
lZ a . Expressing (|3.67|) in terms of (|3.74p we have 



£_(£) =tan" 1 [j_(£ )e 2e( ^ o) 



(3.77) 



where 5-(£q) = sin 1 ^T p (£o) ~ \/2Xq. As before, we 
define X = (T p ° - Tl%)e Ael l2. We may express G Y y = 
(e 2 /h)K a with K a = tt 2 -D 2 . Integrating (f3~76|) we then 
find 



Gyy(X) 



2-Xe~ F ^l\ 
h 



(3.71) where 



F(X) 



-TT 2 , 



2 A' 



dx 



(tan 1 x) ' 



Thus, for X < 0, \X\ > e and e we find 
G 1/2+£ (X) = \X\e~ F ™\ 



(3.78) 



(3.79) 



(3.80) 



The asymptotic behavior F(X) = X/tt 2 for \X\ < 1 
reproduces (|3.58p when \X\ ^> e. For \X\ > 1 we find 

F(X^oo) = Uog2X- 7 -^. (3.81) 

where £(3) = 1.20 is the Riemann zeta function. This 
gives the asymptotic behavior 



Gl/2+e(X -> -OO) = 



; 14C(3)/7r 2 \ 8= 
21X1 



(3.82) 



(3.73) which is quoted in Table HT1 
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IV. DISCUSSION AND CONCLUSION 

In this paper we have examined several novel proper- 
ties of a point contact in a QSHI. We showed that the 
pinch-off as a function of gate voltage is governed by 
a non trivial quantum phase transition, which leads to 
scaling behavior of the conductance as a function of tem- 
perature and gate voltage characterized by a universal 
scaling function. We computed this scaling function and 
other properties of the critical point in certain solvable 
limits which provide an overall picture of the behavior as 
a function of the Luttingcr liquid parameter g. 

In addition, we showed that the four terminal conduc- 
tance has a simple structure when expressed in terms of 
the natural variables, Gab, and that at the low temper- 
ature fixed points, the leading corrections to the differ- 
ent components of Gab can have different temperature 
dependence. In particular, we showed that the skew con- 
ductance Gxy vanishes as T 7 with 7 > 2. 

Finally, we showed that for strong interactions, g < 
1/2, the stable phase is the time reversal breaking insu- 
lating phase. Transport in that phase occurs via novel 
fractionalized excitations that have clear signatures in 
noise correlations. 

There are a number of problems for future research 
that our work raises. We will divide the discussion into 
experimental and theoretical issues. 

A. Experimental Issues 



150nm — this predicts g ~ 0.8. The critical exponent 
governing the temperature dependence of the pinch-off 
curve is then a g ~ .02. In the CC and II phase the 
conductance vanishes as T s with 8 g = g + g^ 1 — 2 ~ .05. 

The good news is that since g is close to 1 the low tem- 
perature scaling behavior should be accurately described 
by the scaling function (|3.9[) computed in the limit g — > 1. 
The bad news, is that the smallness of a g and S g mean 
that it will be difficult to see much dynamic range in the 
conductance as a function of temperature. Nonetheless, 
it may be possible to observe logarithmic corrections to 
the conductance as a function of temperature, and by 
comparing pinch-off curves at different temperatures it 
may be possible to observe the predicted sharpening of 
the transition as temperature is lowered. 

The skew conductance Gxy is predicted be zero for 
non interacting electrons, and with weak interactions 
vanishes as T 2 . This is a consequence of the unique edge 
state structure of the QSHI, and remains robust when 
the interactions are weak. 

To probe the critical behavior of the pinch-off transi- 
tion, as well as the more exotic strong interaction phases 
it would be desirable to engineer structures with smaller 
g. Perhaps this could be accomplished by modifying ei- 
ther the dielectric environment or the bare Fermi veloc- 
ity of the edge states. Maciejko et alii 2 - have suggested 
that this may be possible using InAs/GaSb/AlSb type-II 
quantum well o 46 ' 47 . 



The QSHI has been observed in transport experiments 
on HgTe/HgCdTe quantum well structures. A crucial 
issue is the value of the interaction parameter g. A sim- 
ple estimate can be developed based on the long range 
Coulomb interaction 4 ^. First consider the limit £ 3> w, 
where w is the quantum well width and £ is the evanes- 
cent decay length of the edge state wavefunction into the 
bulk QSHI. We model the edge state as a two dimen- 
sional charged sheet with a charge density profile pro- 
portional to 6(x) exp(— 2x/£), a distance d above a con- 
ducting ground plane. The long range interaction then 
leads to u 2 — u 4 — (2e 2 /e) log(4e 7 ci/£), where e is the 
dielectric constant and 7 = .577 is Euler's constant. As 
a second model, assume £ <C w, and model the edge state 
as a uniformly charged two dimensional strip of width w 
perpendicular to a ground plane a distance d away. This 
gives U2 — 1*4 = (2e 2 /e) log(2e 3 / 2 (i/w). The intermedi- 
ate regime £ ~ w can be solved numerically, and we find 
that it is accurately described by a simple interpolation 
between the above limits with 4d/(£e~ 7 + 2we~ 3 ^ 2 ) in 
the log. This leads to^ 4 -. 



i + i e 



7r efwp 



log 



7. Id 
i + 0.8w 



-1/2 



Kill' 



(4.1) 



30nm 



For e = 15, hv F = .35eVnm, £ = 2hv F /E t 
(i?g a p is the gap of the bulk QSHI), w = 12nm and d 



B. Theoretical Issues 



Our work points to a number of theoretical problems 
for future study. It would be very interesting if the pow- 
erful framework of conformal field theory can be used to 
analyze the intermediate critical fixed point as well as the 
crossover scaling function. Perhaps the first place to look 
is g = l/v3. Maybe it is possible to take advantage of 
the triangular symmetry of the QBM problem to develop 
a complete description of the critical fixed point, analo- 
gous to the mapping to the 3 channel Kondo problem 25 
and the 3 state Potts model^i that apply in a different 
regime. In the absence of an analytic solution, this prob- 
lem is amenable to a numerical Monte Carlo analysis 
analogous to the calculation of the resonance crossover 
scaling function performed in Ref. Il3 . 

In addition, there are a number of other fixed points 
which we did not analyze in detail in this paper. (Recall 
for g = 1 — e we found seven). It would be of interest 
to develop a more systematic classification of all of the 
fixed points, analogous to the analysis of three coupled 
Luttinger liquids performed by Oshikawa, Chamon and 
Affleck and Hou22^. 
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APPENDIX A: FOUR TERMINAL 
CONDUCTANCE 

The electrical response of the point contact can be 
characterized by a four terminal conductance, 

ii=j2 G iM> ( Al ) 

3 

where U is the current flowing into lead i and Vj is the 
voltage at lead j. In this appendix we will develop a 
convenient representation for Gij. Section 1 shows that 
Gij can be characterized by a 3 x 3 matrix, whose en- 
tries have a clear physical meaning. This representation 
allows constraints due to symmetry to be expressed in 
a simple way, which reduces the number of independent 
parameters characterizing the conductance. Finally, in 
section 3 we show how Gij is related to the conductance 
of the SLL model computed by the Kubo formula. 

1. Conductance matrix 

The 4x4 matrix Gij is constrained by current conser- 
vation to satisfy J2 i Gij — J2j Gij = 0. In the absence of 
any symmetry constraints, there are thus 9 independent 
parameters characterizing Gij. In this section we will 
cast these 9 numbers as a 3 x 3 matrix, in which each 
of the entries has a clear physical meaning. In this rep- 
resentation constraints due to symmetry have a simple 
form. 

Since the four currents satisfy ^ - Ii = 0, they are 
determined by three independent currents, which we de- 
fine as I a = {Ix,Iy,Iz), and satisfy 

Ii = J2 M iAlA, (A2) 

a 

where the 4x3 matrix M^a is 

/ i i a 



V i-i-i/ 

Ix = I\ + Ia is the total current flowing from left to right 
along the Hall bar, whereas ly = I\ + h is the current 
flowing from top to bottom. The third current Iz = 
Ii + I3 is the current flowing in on opposite leads (1 and 
3) and flowing out in leads 2 and 4. Similarly, the voltages 



Vi, which are defined up to an additive constant, define 
three independent voltage differences Vp — (Vx, Vy, Vz), 
with 

V B = Y, M BjVj. (A4) 

3 

Vx biases leads 1 and 4 relative to leads 2 and 3, Vy 
biases leads 1 and 2 relative to leads 3 and 4, and Vz 
biases leads 1 and 3 relative to leads 2 and 4. 

The new currents and voltages are then related by a 
3x3 conductance matrix 

i a = J2 GabVb - ( A5 ) 

P 

The 9 elements of Gab determine the four terminal con- 
ductance matrix, 

Gij=Y,M iA G AB Ml j . (A6) 

AB 

The elements of Gab have a simple physical interpreta- 
tion. Gxx is the "two terminal" conductance measured 
horizontally in Fig. [1] by applying a voltage to leads 1 
and 4 and measuring the current 1\ + J4. Similarly Gyy 
is a two terminal conductance measured vertically. Gzz 
describes a two terminal conductance defined by com- 
bining the opposite leads 1 and 3 together into a single 
lead (and similarly for leads 2 and 4). Gxy is a "skew" 
conductance describing the current I\ + 14 in response to 
voltages applied to leads 1 and 2. The other off diagonal 
conductances can be understood similarly. 



2. Symmetry Constraints 

The form of Gab simplifies considerably in the pres- 
ence of symmetries. 



a. Time Reversal Symmetry 

In the presence of time reversal symmetry the four 
terminal conductance obeys the reciprocity relation 42 , 
Gij = Gji. This implies Gab = Gba- Thus, with time 
reversal symmetry the conductance has 6 independent 
components. 



b. Spin Rotational Symmetry 

When the spin S z is conserved the current of up and 
down spins flowing into the junction must independently 
be conserved. It follows that 

Il.in + -^3,in = -^2,out + A.out 

I2, in + Ii, in = h, out + -<3,out- (A7) 
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Since in the Fermi liquid lead (where the interactions 
have been turned off) we have Ii^ n — (e 2 /h)Vi, this im- 
plies that 



h+h = -h -h= -^-(Vi + V 3 - V 2 - Vi). (A8) 
It then follows that 



G zz = 2e 2 /h 
Gzx = Gzy = 0. 



(A9) 



Thus, which spin conservation the conductance is charac- 
terized by 3 components: the two terminal conductances 
Gxx, Gyy and the skew conductance Gxy- 

The quantization of Gzz and vanishing of Gzb are 
therefore a diagnostic for the conservation of spin. 
Though spin orbit terms violating S z conservation are 
generically present, we will argue that at the low energy 
fixed points of physical interest the conservation of spin 
is restored. 



c. Mirror Symmetry 

If the junction has a mirror symmetry under inter- 
changing leads (1,2) <-> (3,4) or (1,4) <-> (2,3), it follows 
that 



Gxy = 0. 



(A10) 



Though mirror symmetry is not generically present in a 
point contact we will argue that that symmetry is re- 
stored in the low energy fixed points of interest. More- 
over, the crossover between the critical fixed point and 
the stable fixed point described by (jl.ip is also along a 
line with mirror symmetry. Thus the crossover conduc- 
tance is characterized by two parameters, Gxx and Gyy, 
which are simply the two terminal conductances. 



d. Critical conductance 

At the transition, where the point contact is just being 
pinched off the two terminal conductances must be equal, 



G 



xx 



G 



YY 



G* 



(All) 



In addition, we will argue that this fixed point also has 
spin rotational symmetry and mirror symmetry. Thus, 
the critical four terminal conductance depends on a 
single parameter G* . 



3. Relation to Kubo conductance 

In this section we relate the conductance matrix Gab 
to the conductances of the SLL model, which can be com- 
puted with the Kubo formula. There are two issues to 



be addressed. First is to translate Gab into the spin and 
charge conductances of the SLL model. Second, we must 
relate the physical conductance measured with leads to 
the conductance computed with the Kubo formula. The 
Kubo conductance describes the response of an infinite 
Luttinger liquid, where the limit L — > oo is taken before 
u> — > 0. This does not take into account the contact re- 
sistance between the Luttinger liquid and the electron 
reservoir where the voltage is defined. An appropriate 
model to account for this is to consider a ID model for 
the leads in which the Luttinger parameter g = 1 for 



x > 



'35,36 



In this section we assume time reversal symmetry and 
that spin is conserved. In this case we may define the 
charge and spin currents in the Fermi liquid leads (x > L) 
to be, 

Ip = Il.'m + I^.in ~ ^l,out — ^4, out 

la = h,in — ^4, in + -fl.out — ^4, out- (A12) 

Similarly, define charge and spin voltages 

V p = (V 1 +V i -V 2 -V 3 )/2 

Vo = {Vi - V 4 + V 2 - V 3 )/2. (A13) 

These are related by the conductance matrix. 

I a = G aP V p , (A14) 



where a, f3 — p, a. By comparing (|A5|) and (|A14[) it is 
clear that 



Gxx 
Gyy 
Gxy 



G pp 

2e 2 /h-G 

G pa — G(jp. 



(A15) 



G a /3 can be computed using the Kubo formula using 
the model in which the interactions are turned off for 
x > L. It is useful, however, to relate this to the Kubo 
conductance G^g of an infinite Luttinger liquid. This 
can be done by relating the voltage V a=P ,a of the Fermi 
liquid leads with g p = g a = 1 to the voltage V a of the 
incoming chiral modes of the Luttinger liquid with g p = g 
and g a — 1/g. By matching the boundary conditions at 
x = L this contact resistance has the form 



with 



V a -V a = R c a0 Ifs 
h g a - 1 



e 2 2g a 



Saf: 



(A16) 



(A17) 



The Kubo formula with infinite leads relates I a — 
G^pVp. Eliminating V a from (|A16|) and (|A17j) gives the 
matrix relation 39 



G 



a/3 



(I - R C G K ) G 



K 



a/3 



(A18) 
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When there is mirror symmetry, so that Gxy = l^po = 
0, the conductance matrix is diagonal, so that (|A18p sim- 
plifies. In that case we find 



Gyy = 2— 
h 



G 



K 

go 



1 R*ooG^na 



(A19) 



Gxx 
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